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Two systematic sources of error in most current microscopic evaluations of fission- 
barrier heights are studied. They are concerned with an approximate treatment of the 
Coulomb exchange terms (known as the Slater approximation) in the self-consistent 
mean fields and the projection on good parity states (e.g., of positive parity for the 
spontaneous fission of an even-even nucleus) of left-right reflection asymmetric intrinsic 
solutions (e.g., around the second barrier). Approximate or unprojected solutions are 
shown to lead each to an underestimation of the barrier heights by a few hundred keV. 



1. Introduction 

In most microscopic calculations of nuclear binding energies using effective nucleon- 
nucleon interaction or their Energy Density Functional (EDF) avatars, data are 
reproduced in the best cases with an accuracy which hardly goes below a couple 
of MeV. When evaluating relative energies, one may hope for a partial cancella- 
tion of such errors but this is by far not granted in general. A particular case of 
choice of such relative energies is constituted by the so-called fission barrier heights. 
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Being a rough one-dimensional description of the very complex many-body fission 
process, their "experimental" values result in most cases from a model dependent 
-or partial- analysis of some data. Nevertheless, they provide undoubtedly a much 
wanted energy scale on which microscopic calculations adjust their phenomenolog- 
ical parameters (e.g. this has been in particular the case for the fit of the SkM* 
Skyrme interaction W). As an example of the accuracy of these barrier calculations 
we may quote the results of Ref. ^ where an average accuracy of about 1 MeV 
for the first and second barriers of actinide nuclei (as well as fission isomer excita- 
tion energies) is found upon making a very rough estimate of the rotational energy 
correction. 

Such a range of error is very significant in physical terms. It is generally esti- 
mated that a variation of 1 MeV roughly corresponds to a change of spontaneous 
fission half lives by 4 orders of magnitude. It is therefore of paramount importance 
to pinpoint systematical errors inherent to a given theoretical (phenomenological) 
approach which could lead to wrong estimates of fission properties or conversely, 
when making a fit of the interaction taking care of fission barriers, which could have 
adverse effects on theoretical estimates of other quantities. 

Without aiming at any exhaustibility in chasing such errors, two examples will 
be studied here: i) the approximate treatment of the Coulomb exchange terms 
known as the Slater approximation in the self-consistent mean fields, ii) the pro- 
jection on good parity states (e.g. of positive parity for the spontaneous fission 
of an even-even nucleus) of left-right reflection asymmetric intrinsic solutions (e.g. 
around the second barrier). 

2. On the effect of the Slater approximation of Coulomb exchange 
terms 

Most microscopic calculations use the infinite matter Pauli correlation function 
(within a local density approximation) to compute the Coulomb exchange expecta- 
tion value in a Slater determinant producing an EDF piece involving only the local 
proton density 

«r e xc h . = ~\e 2 (~) 5 / d\ P y\r) . (1) 
The above yields, within a variational approach, the following mean field piece 

Cr,(r>-(J)'eV 3 «- ( 2 ) 

This approach is usually refered to as the Slater approximation It has been 
tested many years ago ^ by incorporating an exact Coulomb exchange treatment in 
Hartree-Fock calculations of the ground state (spherical or deformed) solutions of 
8 light nuclei (from 16 to 56 Ni). In this early work, the Skyrme SIII interaction^ 
had been used. 
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Such a study has been recently revisited by J. Skalski f or the ground-state 
Hartree-Fock solutions of 9 spherical subshell closed nuclei, from nuclei as light as 
16 to superheavies such as 310 126. The SkP^and SkM* Skyrme interactions have 
been used. 

From these two studies, one has concluded that the errors found (namely rela- 
tive errors, consistently in all this Section) for the Coulomb exchange energies, are 
relatively small, amazingly interaction-independent (at least for the three interac- 
tions under consideration, namely SIII, SkM* and SkP) and decreasing with the 
total particle number A (as expected in view of the infinite matter origin of the 
Slater approximation). 

Even though no quantitative account nor general argument, has been given 
there, as far as fission barrier heights are concerned, one should mention the cal- 
culations of Ref. The authors have included in a HFB framework full Coulomb 
exchange terms in various nuclear systems , including the fission barriers of 254 No. 

In the following we have performed self-consistent (Hartree-Fock or Hartree- 
Fock-plus-BCS calculations) using only the SIII Skyrme effective interaction (with 
a seniority interaction for the BCS part). The Coulomb energy and mean-field 
contributions have been calculated as proposed in Ref. ^ by evaluating the Coulomb 
matrix elements upon using a gaussian integral representation of a Yukawa-type 
two-body interaction as 

= Jl fe-^/V-^lV^ . (3) 

|ri -r 2 \ V 7T J <t 2 

First, we have extended our study of the above spherical subshell closed nuclei 
beyond the Hartree-Fock approximation in use and allowed for pairing correlations a 
la BCS with a seniority force. For two nuclei ( 90 Zr and 298 114) where pair correlated 
solutions have been found to occur, the error has been dramatically reduced (by 
much more than half as seen in Fig. 1). This appears thus as a hint that the error 
found so far is due to the particular filling of those nuclei which had been considered. 

We have then calculated the error for six isotonic series (N = 28, 64, 78, 106, 126, 
146) of nuclei constrained to be spherically symmetric (allowing pairing correlations 
to be included). As seen in Fig. 1 the trend of the error is totally unambiguous. 
Away from proton magic numbers (Z = 20, 28, 50, 82) the quality of the Slater 
approximation improves very rapidly. To confirm the role of the proton single- 
particle level scheme, we have performed similar calculations for three short isotopic 
series (Z — 72, 82, 92). The error is found, see Fig. 1, to be remarkably stable within 
a given isotopic series. 

Then, we have explored changes in the proton single-particle level density pp 
at the Fermi surface, not upon changing the proton number but instead the axial 
quadrupole deformation for the same nucleus. As examplified in Fig. 2 reporting 
the results of calculations for 5 nuclei ( 24 Mg, 48 Cr, 106 Mo, 178 Hf, 238 U) the error is 
always-and often very significantly-larger for their deformed ground-state solutions 
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Fig. 1. Relative error Ai?cE associated with open proton shell nuclei, where Ece is the expectation 
value of the exchange part of the Coulomb interaction. The closed proton shell nuclei are reported 
as open circles linked by the solid line. Some isotonic series are represented in dash lines, whereas 
three isotopic series corresponding to Hf, Pb, and U isotopes are plotted with black triangles. For 
the 90 Zr and 298 H4 nuclei, the plotted solutions correspond to BCS solutions. 



than for the solutions constrained to have a vanishing quadrupolc moment. In the 
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Fig. 2. Relative errors Ai?cE in the ground-state solution of deformed nuclei and in their spheri- 
cally constrained solution. The solid line, shown for comparison, corresponds to the closed proton 
shell nuclei appearing in Fig. 1. 
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case of 238 U, we have also reported the errors corresponding to the top of the first 
barrier and the fission isomeric state. These errors lie between those obtained for 
the spherical and deformed solutions. To confirm these findings we have compared 
the error made on the Coulomb exchange energies (middle panel of Fig. 3) with 
the proton BCS pairing gap, taken as an index of the proton single-particle level 
density at the Fermi surface (lower panel of Fig. 3) along the beginning of the 
calculated fission barrier of the 70 Se nucleus imposing a vanishing axial octupole 
moment (upper panel of Fig. 3). As seen from Fig. 3 the error and the gap are very 




Q20 [ fm2 l 



Fig. 3. (Top panel) Deformation-energy curve of the 70 Se as a function of the mass quadrupolc 
moment Q20 in fm 2 2. (Middle panel) Relative error AEce as a function of Q20- (Bottom panel) 
Variation with Q20 of proton pairing gaps calculated with the exact treatment of Coulomb ex- 
change terms in solid line (with full circles) and using the Slater approximation in dashed line 
(with open circles). 

well anti-correlated functions of the deformation. 

We can therefore conclude that the quality of the Slater approximation will be 
less good for low values of than for densities corresponding to mid-shell situations. 

Since it is well known that systematically p will be higher at the top of the 
fission barrier than at the ground state deformation (for a direct evidence of that 
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fact in the context of fission-barrier calculations within the Hartree-Fock-plus-BCS 
framework, see Ref. \ we can predict that the correction to be made to account 
for the error brought in by the Slater approximation will be lower in the former case 
than in the latter. As a result, using the Slater approximation will systematically 
lead to an underestimation of the fission-barrier he igh ts. In the calculations already 
discussed (see Fig. 2) of the 238 U fission barrier one should add beyond the 
Slater approximation values about 310 keV to the first fission barrier height and 
about 280 keV to the fission isomeric energy, quantities which are not without 
consequences on spontaneous fission halflives as already mentioned. Even though 
it is hard to get more than a tendency on the l.h.s part of the Figure 12 of Ref. ^ 
their results (for 254 No) are consistent with our above conclusions. 

Before closing this Section, two facts deserve a short notice. 

In the above discussion of the impact of ~p^, on the quality of the Slater approx- 
imation, we have only considered the errors on the Coulomb exchange energies. Of 
course, the quantities whose errors should have been compared (and actually which 
have been discussed just above in the 238 U fission study) are total energies (the 
sum of kinetic, nuc lear , direct and exchange Coulomb energies). However, it has 
been demonstrated that, in most cases, the errors on such total energies and 
Coulomb exchange energies are equal up to some tens of keV. 

As found out already in the results of Ref. El and emphasized in Ref. it has 
been systematically found that all hole states, in Hartree-Fock parlance, are pushed 
up by the Slater approximation from their location in exact calculations, while 
conversely, all particle states are pushed down. As a result the proton single particle 
level density at the Fermi surface is unduly enhanced by the Slater approximation. 
The effect of such an error on proton particle-hole energies as far as correlation 
properties are concerned, remain to be specifically studied. Indeed, if the residual 
interaction is to be considered as a perturbation, it is easy to see that such a trend 
could affect significantly the configuration mixing. 



3. On the effect of a projection on good parity states 

For a very long time, it is known that the second fission barriers of actinides 
are unstable with respect to the left-right asymmetry. It has been then recognized 
that this phenomenon is driven by the heavy fragment shell effects explaining thus 
already at this early stage of the fission process, the asymmetric pattern of the 
fragment mass yields (at low compound-nucleus energy) which had been observed 
in this region long before. 

In most microscopic or macro- microscopic calculations in the U, Pu region, upon 
increasing the elongation after the fission isomeric state, the intrinsic equilibrium 
solution becomes unstable with respect to the left-right reflection symmetry and 
acquires rapidly a larger and larger octupole deformation which stabilizes at a value 
corresponding, as we just have noted, to the most probable fragmentation. This fact, 
as we will see, plays an important role for our purpose here. 
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Yet, even though the intrinsic parity may be broken for some microscopic so- 
lutions under consideration, the parity of the physical solution must be conserved 
during the fission process. For instance, if one describes the spontaneous fission of 
an even-even nucleus as 240 Pu, one should evaluate the fission barrier obtained upon 
imposing the positive parity to the solutions. This may be obtained by projecting 
intrinsic solutions on the desired parity. This corresponds thus to a projection after 
variation, amounting to a mixing of two configurations and |$) = IT | M/) , where 
n is the parity operator. 

Recently such a projection of correlated microscopic solutions a la HTDA 
(Highly Truncated Diagonalization Approach) ^ has been extensively studied as 
the PhD thesis work of one of the authors We will thus simply outline here, 
the major points of this approach and discuss the results which have been obtained 
when applying it to the second fission positive-parity barrier of Pu. 

The HTDA calculations are designed to produce realistic correlated wave func- 
tions in an approach which conserve explicitly the particle number and do not 
violate the Pauli principle. They consist in performing intrinsic shell-model like cal- 
culations using single-particle states deduced from a mean-field potential U. This 
field is a priori arbitrary. Yet, its realistic character, given a state-of-the-art micro- 
scopic effective Hamiltonian H , is essential to limit the size of our n-particle n-hole 
many-body basis. The mean field U will be taken here as the one-body reduction 
of the nucleon-nucleon effective interaction V included in H for a corresponding 
approximate (e.g., a la BCS) microscopic correlated wave function. 

The HTDA Hamiltonian is written as 

H = K + V = H MF + Vrcs (4) 

where K is the kinetic energy while the mean field (one body) Hamiltonian Hmf 
is defined by 

Hmf = K + U-C (5) 
and the residual interaction Vi-es by 

V Iea = V-U + C (6) 

The particle- hole vacuum |$q) is defined as an eigensolution of Hmf- Here we 
have chosen this vacuum as the ground state of Hmf- The constant C appearing 
in the above definitions is chosen in such a way that the expectation value of the 
residual interaction is vanishing for |<&o)- 

In the present calculations, we have restricted the particle and hole single- 
particle states to lie within 6 Mev below and above the Fermi energy (defined as half 
the sum of the energies of the last occupied and first unoccupied state in |3>o))- The 
many body basis is comprised, beyond the vacuum, of only one-pair-transfer states 
(pairs meaning here Cooper pairs of time-reversal conjugate single-particle states). 
We have chosen the Skyrme SkM* for V. As usually done in practical HTDA cal- 
culations, we have replaced V by a density-independent zero-range interaction to 
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define the residual interaction. The correlated wave function \^f) is obtained by a 
diagonalization of the Hamiltonian H through a Lanczos al gori thm process. 

Through a judicious choice of the mean field U (see Ref. EH) the expression for 
the energy of the projection of the state |\E') on a state of good parity p = ±1 may 
be simplified as 

The non-diagonal overlaps are calculated using methods due to Lowdin (The 
use of the heavy Balian-BrezinEIl generalized Wick theorem is not necessary in our 
case since we deal only with Slater determinants). 

The results of our projection (after variation) calculations for the fission barrier 
of 240 Pu for a positive parity, are summarized in Figure 4. In the lower panel, it 
appears that the projection has the following effects: i) it yields a small stable 
octupolc deformation around the fission isomeric state, defined by an elongation 
Q20 , f° r an elongation much smaller than the one (Q20 ) where an instability of the 

HTDA unprojected solution is observed, ii) before reaching the second (asymmetric) 

(2) 

fission barrier at an elongation Q 2Q , the projection does not yield any significant 
effect on the energies. 

This can be explained in the following way which is illustrated in the upper panel 
of Fig. 4. Much before Q 2 o , a ^ Q20 = 64 b for example, the unprojected solution is 
stable with respect to the axial Q30 mode with a stiffness large enough to only allow 
for the creation, by the projection on positive parity, of a shoulder on the E(Q 30 ) 
curve at a fixed elongation. Around the superdeformed solution Q 2 °q ~ 85 b the 
stiffness parameter is no longer sufficiently large to prevent the projection from 
creating a stable equilibrium deformation (such a weak stiffness had already been 
found in Hartree-Fock-Bogoliubov calculations in the super-deformation region of 
Hg and Pb isotopes G^j). 

At a quadrupole deformation around Q20 ~ 110 b, the non-diagonal overlaps (for 
the identity and Hamiltonian operators) become small enough so that the projected 
energies are not significantly different from the unprojected ones, a property which 
will be, of course, all the more verified than the octupole deformation will increase. 
For even larger elongations, e.g., Q20 ~ 130 b, the equilibrium octupole deformation 
parameter will get larger to stabilize, as we have recalled, at a value typical of the 
most probable asymmetric fission fragmentation and thus yielding a minimum in 
the E{Qso) curve exactly at the same place than in the unprojected case. 

As a result the height of the second fission barrier when projected on the positive 
parity will be enhanced from its unprojected value. This effect is only due to the 
relatively soft character of the octupole deformation energy curve near the fission 
isomeric state and not to the projection at the top of the second barrier where the 
octupole deformation is too large to yield any projection effect. 

It is our contention that these two properties (for the isomeric state and at the 
top of the second fission barrier) are quite general in the actinide region. Therefore, 



4, 2012 0:46 WSPC/INSTRUCTION FILE Quentincorrected 





Fig. 4. Deformation energy curves for the 240 Pu nucleus in the fission isomeric state and second- 
barrier region. Upper panel: energy curves as functions of the axial octupole moments for fixed 
values of the axial quadrupole moment. Lower panel: resulting fission barriers. 



we deem that the underestimation of such fission-barrier heights is probably a 
systematic effect. In our current calculations, which are somewhat approximate, 
this enhancement of the fission-barrier height is of about 350 keV. 

4. Conclusions 

We have discussed here two systematic effects leading to the conclusion that usual 
microscopic calculations of actinide fission-barrier heights could be underestimated 
by a few hundred keV for each. 

One effect is related to the correlation existing between the quality of the so- 
called Slater approximation for the Coulomb exchange energy and the value of the 
proton single-particle level density. En passant, we note that this defect (being 
related to a purely quantal property) could not be possibly cured by a fit of EDF 
parameters. Moreover it induces a systematic lowering of the proton particle-hole 
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energy which might play a substantial role in the correlations built upon using such 
approximate single-particle states. 

The second effect has been established from parity-restoration calculations from 
HTDA correlated solutions. Beyond the systematic conclusion they seem to indi- 
cate, they correspond to the first configuration-mixing calculations of such HTDA 
states where the advantage of their Slater determinantal character has been fully 
exploited. 

References 

1. J. Bartel, P. Quentin, M. Brack, C. Guet and H.-B. Hakansson, Nucl. Phys. A 386 
(1982) 79. 

2. L. Bonneau, P. Quentin and D. Samsoen, Eur. Phys. J. A 21 (2004) 391. 

3. J.C. Slater, Phys. Rev. 81 (1951) 385. 

4. C. Titin-Schnaider and P. Quentin, Phys. Lett. B 49 (1974) 397. 

5. M. Beiner, H. Flocard, Nguyen Van Giai and P. Quentin, Nucl. Phys. A 238 (1975) 
29. 

6. J. Dobaczewski, H. Flocard and J. Treiner, Nucl. Phys. A 422 (1984) 103. 

7. M. Anguiano, J.L. Egido ansd L.M. Robledo, Nucl. Phys. A 683 (2001) 227. 

8. P. Quentin, J. Phys. (France) 33 (1972) 457. 

9. M. Brack and P. Quentin, Nucl. Phys. A 361 (1981) 35. 

10. J. Le Bloas, Meng-Hock Koh, P. Quentin, L. Bonneau and J. I. A. Ithnin, Phys. Rev. 
C84 (2011) 014310. 

11. J. Skalski, Phys. Rev. (7 63 (2001) 024312. 

12. P. Moller and S.G. Nilsson, Phys. Lett. B 31 (1970) 283. 

13. N. Pillet, P. Quentin and J. Libert, Nucl. Phys. A 697 (2002) 141. 

14. T.V. Nhan Hao, Ph.D. Thesis, Bordeaux 1 Univ., Talence, France (2010). 

15. P.O. Lowdin, Phys. Rev. 97 (1955) 1474; 1490. 

16. R. Balian and E. Brezin, Nuovo Cimento B 64 (1969) 37. 

17. J. Skalski, P.-H. Heenen, P. Bonche, H. Flocard and J. Meyer, Nucl. Phys. A 551 
(1993) 109. 



